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ABSTfiACT 


The  time-averaged  power-density  distribution  below 
the  ocean  surface  due  to  incident  laser  radiation  is 
examined  by  means  of  computer  simulation  of  the 
geometrical  optics  involved  with  the  air/sea  interface 
and  subseguent  ocean  penetration  by  the  laser  beam. 
The  effects  over  the  entire  spectra  of  incidence 
angles,  wind  velocities,  wind  directions,  beam  spot 
sizes  and  depths  of  penetration  are  analyzed. 
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I.  I NT BOD OCT I ON 


The  tine-averaged  power-density  distribution  below  the 
ocean  surface  due  to  incident  laser  radiation  is  examined 
over  the  entire  spectra  of  incidence  angles,  wind  velocities 
and  wind  directions  by  means  of  geometrical  optics. 

Current  interest  in  using  laser  systems  for 
communications  through  the  air-sea  interface  and  for 
detecting  submerged  objects  necessitates  the  development  of 
prediction  methods  for  determining  the  subsurface 
power-distribution  of  radiation  from  a laser  source  above 


ocean. 


mathematical 


model 


optical-communications  application  done -by  Karp  in  Bef.  1 is 
developed  in  terms  of  a radiance  function  related  to  the 
mutual  coherence  function  with  accurate  results  for 
incidence  angles  out  to  45  degrees.  A more  general  model 
developed  by  Swennen  in  Befs.  2-4  examines  the  power-density 
distribution  below  the  ocean  surface  through  a rigorous 
assessment  cf  the  ocean  surface  geometry  allowing  for 
theoretically  accurate  results  over  the  entire  incidence 
angle  spectrum. 


The  computer  code  developed  for  this  analysis  of  the 
power-density  distribution  below  the  ocean  surface  is  an 
expansion  of  Swennen*s  theory[Befs.  2-4]. 


6 


II.  DESCRIPTION  OF  PHYSICAL  PARAMETERS 


A.  OCEAN  SURFACE 


A geometrical  representation  of  the  ocean  surface  in 
terms  of  local  slopes  developed  by  Cor  and  MunJt[Ref.  5]  is 
the  basis  for  the  analysis  of  the  power  probability 
distribution  below  the  ocean,  as  it  is  in  Refs.  2-4.  A 
review  of  this  representation  follows. 

The  center  of  symmetry  of  the  incident  laser  beam  cn  the 
ocean  surface  facet  is  the  center  of  a right-  handed 
cartesian  coordinate  system,  point  "O"  in  Fig.  1.  The 
z-axis  points  vertically  upward.  The  y-axis  is  in  the 
horizontal  plane,  colinear  with  the  projection  of  the  center 
of  the  incident  laser  beam  onto  that  plane  and  pointing  in 
the  direction  of  the  laser  source. 

The  slopes  of  the  ocean  surface  facet  are  defined  by 
angular  parameters  Alpha  and  Beta.  The  angle  Beta  is  the 
angle  between  ^.the  line  of  steepest  ascent  of  the  facet  and 
the  x-y  plane;  Alpha  is  the  angle  between  the  projection  of 
the  line  cf  steepest  ascent  onto  the  x-y  plane  and  the 
y-axis. 


A second  cartesian  coordinate  system  at  the  surface  is 
used,  as  explained  in  Ref.  2,  to  simplify  the  computation  of 
the  slope  probability  function.  This  coordinate  system  is 
designated  by  primes  with  the  y'-axis  pointing  towards  the 
wind  source  and  rotated  an  angle  Chi  from  the  y-axis.  The 
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z'-axis  coincides  with  the  z-axis.  The  slope  parameters  are 
also  primed  with  their  transformations  given  by: 


Eeta'  = Beta 
Alpha*  = Alpha  - Chi 


I ; 


Pi 


i 


■ ! 
► 1 


Belcw  the  ocean  surface  a depth  Z at  the  point  of 
observation.  O',  is  centered  a translational  transformation 
of  the  surface  coordinate  system.  Angles  Mu  and  Nu  define 
the  laser  beam  refracted  to  0 1 by  any  surface  facet.  The 
first  angle.  Mu,  is  the  angle  between  the  refracted  ray  and 
the  z-axis;  the  second,  Nu,  is  the  angle  between  the  y-axis 
and  the  vertical  plane  of  the  refracted  ray. 

B.  LASER  BEAM 


A cylindrical  laser  beam  of  radius  a is  incident  on  the 
ocean  surface  wave  facet  at  an  angle  Psi  measured  from  the 
z-axis.  The  intersection  of  the  beam  and  the  facet  is  a 
flat  elliptical  surface  as  illustrated  in  Fig.  2.  The 
semi-minor  axis  of  the  ellipse  is  equal  to  the  beam  radius 
and  the  semi-major  axis,  b,  is  defined  by: 

^ \ 

b = a * sec  Psi 


The  analytical  development,  as  in  Refs.  2-4,  begins  at 
the  ocean  surface  facet.  The  laser  source,  aiming, 
atmospheric  propagation,  etc.  are  extraneous  to  this 
analysis.  Pcwer  density  values  at  the  surface  are  not 
quoted  because  the  analysis  deals  with  the  ratio  of  the 
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power  density  at  the  depth  of  interest,  P , and  the  power 

d 

density  at  the  surface,  P 

ds 

C.  MAXIMUM  ANGLE  OF  INCIDENCE 

The  angle  of  incidence  with  respect  to  the  z-axis,  Psi, 
of  the  laser  beam  is  analyzed  from  zero  degrees  to  a maximum 
of  84  degrees.  At  greater  angles  the  beam  will  not  be  able 
to  strike  all  the  wave  facets  within  a wavelength  cf  the 
wave  due  to  the  wave  height.  There  is  also  the  possibility, 
as  illustrated  in  Fig.  3,  that  a ray  incident  at  Psi  > 84 
degrees  could  pass  through  the  wave,  re-enter  the  atmosphere 
and  then  re-enter  the  ocean;  or,  depending  on  the  local 
slopes,  the  ray  could  undergo  a total  internal  reflection  if 
the  Critical  Angle  were  exceeded.  The  mathematical  model  ‘of 
Swennenf Refs.  2-4]  and  this  thesis  do  not  address  this 
situation. 

This  limitation  on  the  maximum  angle  of  incidence, 

Psi  , is  derived  from  Fig.  4 and  defined  by: 
max 

-1 

Psi  = tan  (L/2H) 
max 


The  significant  wave  height,  H,  and  the  average  period,  T, 
of  the  wave  are  given  by  equations  (1)  and  (2)  as  defined  by 
PiersonfBef.  7].  The  average  wavelength  of  the  wave  is 
given  by  equation  (3)  as  defined  by  HillfRef.  8]. 


■wwMwammiiwwwm^iMWr  » - - 


a 


1 = 0.81  * (2H/g) 
2 

1-91/2 


where 

W = wind  velocity 
g = gravitational  constant 

For  wind  velocities  of  1 to  14  m/sec  the  maximum 
of  incidence  is  84.2  degrees  (±0.4  degrees)  for  all 
velocities;  therefore,  Psi  is  limited  to  84  degrees 


(2) 

(3) 


angle 

wind 


Ill 


iOjill-DENSITY  PROABILITY  INTEGRAL 


A.  GENERAL  EQUATION 


The  power-density  probability  distribution  (P  ) at  the 

d 

point  of  observation  below  the  ocean  surface  (O')  due  to  the 

entire  beam  is  given  by  the  integral  of  dP^  over  all  of  the 

contributing  surface  facets.  As  a function  of  the  angular 

parameters  (Mu,Nu)  the  integral  as  defined  in  Refs.  2-4  is 
given  by: 


where 


P 

d 


(Mu, Nu) 


dMu  dNu 


F (Mu, Nu)  = P *t*DTF*cos (HI) *sec  (US) 
ds 


*[  cos  (Mu)  *cos  (Mu^)  +sin  (Mu)  *sin  (Mu^)  *cos  J Nu-Nu^  | ] 

2 -1 

*P  (Z* , 2') *tan (Beta) *sec  (Beta) *J 
x y 

and 


E = power  density  at  the  surface  facet 
ds 

T = Fresnel's  transmittance  coefficient 
ETF  = diffuse  transmittance  function 
HI  = angle  of  incidence 
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HE  = angle  of  refraction 

Mu^  = first  quadrant  angle  the  refracted  center 


ray  of  the  beam  makes  with  the  positive 
z-axis 

= positive  angle  the  projection  of  the 


refracted  center  ray  of  the  beam  onto  the 

xy-plane  makes  with  the  y-axis  measured 

from  the  latter  to  the  former 

F(Z',Z')  = time-average  slope  distribution 
x y 

function 


The  development  of  equation  (4)  is  described  in  great 
detail  in  Refs.  2-4  and  will  not  be  repeated  here;  however, 
for  clarity  the  major  components  of  equation  (4)  will  be 
discussed. 


i 

B.  FRESNEL'S  TRANSMITTANCE  FUNCTION 


Fresnel's  transmittance  coefficient  (T)  is  defined  as 
the  intensity  ratio  of  the  transmitted  to  the  incident  beam 
energy. 


The  collimated  incident  laser  beam  is  composed  of  two 
polarization  components,  a transverse  magnetic  or  parallel 
polarization  (T^)  and  a transverse  electric'  or  normal 


polarization 


(T  ) 

1 


Resolving  the  refracted  beam  into  its 


components  (T  ,T  ) is  a complex  task  due  to  rotation  of  the 
l // 

plane  of  incidence  about  all  three  coordinate  axis  over  the 
range  of  integration.  In  order  to  deal  with  this  complexity 
the  beam  is  assumed  to  be  rendered  diffuse  and  unpolarized 
below  the  ocean  surface  due  to  scattering,  propagation 
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direction  variations  and  phase  changes.  This  assumption 

allows  for  the  averaging  of  T and  T to  obtain  the  value 

1 // 


of  I. 


T = (T  +T  )/2 
1 // 


The  equations  for  T and  T used  by  Swennen[  fief s . 

1 // 


are: 


2-4] 


2 

T = [ sin  (MI)  *sin  (Wfi)  /sin  (Wl+Wfi)  ] (5) 

1 

2 

•I  = C sin  (WI)  *sin  (HH) /sin  (WI+WH)  *cos  (HI-WR)  ] (6) 

Equations  (5)  and  (6)  are  in  disagreement  with  those 
developed  by  Ecrnf fief . 9]  and  Fowles[Ref.  10],  given  by: 

2 

T = i 2*cos  ( WI)  *sin  (W3)/sin  (WI  + HR)  ] (7) 

1 

2 

I = [ 'os(WI)  *sin(«R)//sin  (HI  + WR)  *cos(Wl-WR)  ] (8) 

Equations  (5)  and  (6)  will  give  erroneously  high  values 
for  the  Fresnel  transmittance  function,  especially  at  higher 
angles  of  incidence.  Because  of  this,  equations  (7)  and  (8) 
are  used  in  this  modeling. 

The  geometry  of  the  ray  refraction  through  the  ocean 
surface  is  illustrated  in  Fig.  5. 


i 


13 


R 


C.  DIFFUSE  TRANSMITTANCE  FUNCTION 


The  diffuse  transmittance  function  (DTF)  accounts  for 
beam  attenuation  below  the  ocean  surface  due  to  scattering 
and  absorbticn. 

The  six-constant  DTF  developed  by  Duntley  in  Ref.  1 1 is 
simplified  tc  a function  of  two  constants;  the  backward 
scattering  coefficient  (BS)  and  the  total  absorbtion 
coefficient  (AC) . This  simplification  is  valid  because  of 
the  random  scattering  particle  orientation  in  the  ocean. 
The  DTF  is  defined  by: 


DTF  = K/[  (AC+BS)  *sinh  (K*Z) +K*cosh  (K*Z)  ] 
where 

V2 

K = [ AC*  (AC+2BS)  ] 

Z = depth  of  the  point  of  observation 

In  current  terminolgy  this  method  of  accounting  for 
attenuation  is  a zero-angle  forward  scattering  technique. 


Typical  values  for  BS  and  AC  were  used  from 
measurements  by  Tyler[Ref.  12]  and  Duntley[Ref.  13].  These 
represent  an  average  between  cool  and  warm  ocean  waters  in 
the  blue-green  spectrum  (4800  angstrom). 


BS  = 0.065  percent 


AC  = 0.044  m 
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D.  SLOPE  PECEABILITY  FUNCTION 


The  time-average  slope  probability  function,  P(Z',Z'), 

x y 

is  a statistical  distribution  of  the  ocean  surface  slopes 

developed  by  Cox  and  Munk[Ref.  5].  The  distribution  derived 
from  the  surface  geometry  is  Gaussian,  which  is  then  altered 
by  a Gram-Charlier  series  in  order  to  account  for  slope 
skewness  and  peakedness  caused  by  the  wind. 

P (Z',Z*)  = f (Xi,Eta,M) 

x y 

where 

Xi  = standardized  crosswind  slope  component 
Eta  = standardized  upwind  slope  component 
H = wind  velocity 

The  development  and  the  equation  for  the  slope 
probability  function  are  covered  in  detail  in  Refs.  2,4,  and 
5.  A clean  ocean  surface  is  assumed  and  the  limits  of 
applicability  placed  on  the  functional  parameters  are 
adhered  to,  namely: 

I Xi | < 2.5 
| Eta  j < 2.5 
N < 14  m/sec 


E.  JACO  El AN 


The  power-density  probability  integral,  equation  (4) , 
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was  first  developed  over  the  slope  components  (Z',Z').  It 

x y 

was  then  transformed  to  a function  of  the  angular  parameters 

(Alpha, Beta)  and  then  by  the  Jacobian  (J)  to  a function  of 
(Mu,Nu)  . 

■ fK 

z;z; 

= |"Jf  (Alpha, Beta)  dAlpha  dBeta 
= J"  Jf  (Mu,  Nu)  J 1 dMu  dNu 

where 


J = ^(Mu,Nu)/^(  Alpha,  Beta) 
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IV.  METHODS  OF  SOLUTION 


A.  FAB  ZONE 


Swenaen[ Refs.  2-4]  observed  that  when  the  depth  (Z)  of 
the  point  of  observation  is  large  compared  to  the  beam 
cross-sectional  radius  at  the  surface  (Z/a>100)  the 
integrand  of  equation  (4)  remains  essentially  constant 
during  the  integration.  The  integral  can  then  be 
approximated  by: 


= J[f(Mu, 


Nu)  *J  dMu  dNu 


- F(Mu,Nu)*J  a Mu  a Nu  ; with  Mu=Muq  & Nu=Nuq 


P = p *T*tan (Beta) *sec  *cos(HI) 
d ds 


♦sec  (WR)  (Z  ' # P' ) *»7  ♦aA 

* y 


The  area  aMuaNu  in  the  (Mu,Nu)  plane  is  approximated  by 
an  ellipse  of  area  aA,  whose  development  is  covered  in 
detail  in  Ref.  3. 

This  is  the  Far  Zone  approximation  and  equation  (4) 
becomes: 


I I • | ■ — ■ 


The  basic  Far  Zone  computer  program  developed  in 

FORTRAN  language  and  run  on  the  IBM  360  system  solves  for 

the  ratio  of  the  power  density  at  the  point  of  observation 

to  the  power  density  at  the  surface  facet,  P /P  . The  main 

d ds 

portion  cf  the  program  allows  for  the  entry  of  any 

combination  of  variables  (wind,  depth,  beam  radius,  Psi, 

Chi,  Hu  , backward  scattering  coefficient,  absorbtion 
0 

coefficient)  , and  also  computes  the  diffuse  transmittance 

function,  Fresnsl* s transmittance  function  and  the  slope 
probability  function.  Subroutines  common  to  all  variable 
entries  are  used  to  compute  the  angles  (Alpha,  Beta,  VI, 
VR) , the  Jacobian  and  A. 

It  is  possible  to  investigate  the  power-density 

probability  distribution  over  a wide  range  of  variables  with 
a minimal  cf  storage  requirements  and  computation  time  due 
to  the  use  of  an  IBM  System/360  Source  Library  subroutine, 
HLHSYS,  in  solving  for  the  coupled  (Alpha, Beta)  and  (Mu, Hu) 
angle  pairs  from  the  simultaneous  non-linear  equations: 

cos  (Mu  ) = [cos  (Psi)  ♦ K*cos  (Beta)  ]/n  (9) 

0 

cot  (Nu  ) = cot (Alpha)  - sin  (Psi) /K 
0 

♦sin (AlphA) *sin (Beta)  (10) 


where 


n = relative  index  of  refraction  (1.33) 
K = K (Psi, Alpha, Beta, n) 
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J 

B.  SEAS  ZCNB 


When  the  depth  or  the  point  of  observation  is  not  large 
compared  to  the  beam  cross-sectional  radius  at  the  surface 
(Z/a<100)  the  Par  Zone  approximation  does  not  hold  and  the 
double  integration  of  eguation  (4)  must  be  carried  out. 

Swennenffief.  2]  used  Simpson's  method  of  numerical 
integration  to  approximate  equation  (4) . The  method  was 
subject  to  singularity  points  in  the  integration  limits 
which  increased  the  complexity  and  time  required  for  the 
computation . 

The  double  integration  can  be  more  accurately  and 
relatively  simply  approximated  by  using  the  Gauss  Quadrature 
system  of  solution,  Eefs.  14  and  15.  The  Two-Point 
Gauss-Legendre  Quadrature  method  was  used.  This  method, 
valid  for  polynomials  up  to  degree  three  (equation  (4)  is  of 
degree  two) , consists  of  first  transforming  the  function 
F(Mu,Nu)  into  a function  F(s,t)  whose  interval  is  -1<s<1  and 
-1<t<1  by  letting: 


Mu  = £ (Mu  -Mu  ) *s+Mu  ♦ Mu  ]/2 
U 1 u 1 

(11) 

Nu  = [ (Nu  -Nu  )*t  + Nu  + Nu  1/2 
U 1 u 1 

(12) 

and 


dMu 

= £ (Mu  -Mu  ) /2  Ids 

(13) 

U 

1 

dNu 

* £ (Nu  -Nu  ) /2  ldt 

(14) 

u 

1 

, 

1 

i 
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The  power-density  integral 
>i-Vn 

P^  *Jjp(Mu,Nu)  dMu  dNu 
«l‘Ul 


is  transformed  into 


P = [ (Mu  -Mu  ) (Nu  -Nu  )/4]*  (fp(s,t)  ds  dt  (15) 
d u 1 u 1 J J, 

-1  *i 

The  Two-Point  Gauss-Legendre  Quadrature  approximation 
of  equation  (13)  is  given  by: 

1 1 

P = [ (Mu  -Mu  ) (Mu  -Nu  )/4  ^ B *F(s  ,t  } (16) 

d u 1 u 1 j*oj*o  i 1 j 


vbere 


W - 1.0  - Gauss  weight  factors 
i 


s.,t^  = ± (3)  = Gauss  roots 


The  basic  Near  Zone  computer  program  developed  uses 

this  Gauss-Legendre  Quadrature  method  in  solving  the 

power-density  integral  coupled  with  the  same  subroutines  as 

the  Par  Zone  program.  The  main  program  computes  the 

integration  limits  (Mu  ,Mu  ) and  (Nu  ,Nu  ) after  determining 

u 1 u 1 

whether  the  surface  projection  of  the  point  of  observation 
is  inside  or  outside  the  beam's  horizontal  cross-secticnal 
area.  The  equations  for  the  integration  limits  are  well 
defined  in  Refs.  2-4. 


The  transformations  described  by  equations  (11)- (15) 
are  then  performed  and  the  power-density  is  solved  for  using 
equation  (16).  As  in  the  Far  Zone  program,  the  Near  Zone 
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program  can  take  any  combination  of  variables  as  entries. 


7 . RESULTS 


Results  for  oblique  incidence  angles  (0<Psi<84  degrees) 
were  obtained  for  vertical  plane  cuts  at  Nu^  = 0 degrees 

over  the  entire  range  of  variables;  however,  valid  data  for 

vertical  planar  sections  at  Nu^  * 0 degrees  could  not  be 

obtained  due  to  values  of  the  standardized  crosswind  and 


upwind  slope  components  being  out  of  the  range  of 

applicability  for  the  slope  probability  function.  This 
problem  is  believed  to  be  inherent  in  the  method  of 
mathematical  analysis  of  the  ocean  surface  used.  Results 
for  normal  incidence  (Psi  = 0 degrees)  over  the  entire  range 


of  variables, 
obtained. 


including  Nu  planar  variations,  were 


The  validity  check  carried  out  on  the  slope  probability 
function,  as  discussed  in  Section  III.D,  clearly  identified 
invalid  data.  Another  source  of  invalid  data  occurred  at 
higher  incidence  angles  (Psi)  where  the  simultaneous 
solution  of  ncn-linear  equations  (9)  and  (10)  resulted  in  an 
angle  of  incidence  (WI)  greater  than  90  degrees  while  the 
slope  probability  function  still  indicated  valid.  This 
output  was  also  easily  identified  and  eliminated. 

Computations  with  the  Near  Zone  solution  method  used 
aore  computer  time  and  storage  than  did  the  Far  Zone 
solution  method.  For  a computer  solution  using  one  value 
for  each  variable  entry,  the  Near  Zone  used  36  per  cent  more 
computer  time  and  8 per  cent  more  computer  storage  than  did 
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1 


’ 


. 


< 


: 1 


an  equivalent  Far  Zone  solution.  The  higher  time 
requirement  for  the  Near  Zone  solution  was  due  mainly  to  a 
greater  compiling  time  requirement.  The  time  difference 
between  the  two  solution  methods  decreased  as  the  number  of 
variables  that  were  incremented  during  one  computer  run 
increased. 

A sample  computer  run  with  output  is  presented  in 
Appendix  E. 


A.  PRESENTATION  OF  RESULTS 


Results  are  presented  graphically  to  simplify 
comparisons  cf  the  wide  range  of  variables.  Vertical-planar 
sections  cf  the  power  density  probability  distribution  below 
the  ocean  surface  cut  through  the  center  of  the  beam 
cross-section  at  the  horizontal  ocean  surface  at  an  angle 

Nu^  are  used  for  the  majority  of  the  results  presented. 

This  presentation  is  in  rectangular  coordinates  with 

normalized  power  density  in  decibels  on  the  ordinate  and  the 
angle  Mu^  in  degrees  on  the  abscissa.  The  normalized  power 

2 

density  is  equal  to  10*log  [ (P  *Z  ) /P  ].  This 

10  d ds 

. - . 2 

normalization  with  respect  to  P /z  should  produce  a 

ds 

distribution  function  that  is  independent  of  depth  in  the 

Far  Zone  regime,  according  to  Swennen[Ref.  2-4],  because  the 
multiple  slopes  seen  in  a beam  cross-section  at  the  surface 
tend  to  render  the  radiation  diffuse.  The  power  density  of 
diffuse  radiation  at  large  depths  decreases  as  the  square  of 
the  depth.  When  the  effects  of  scattering  and  absorption 


I 
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are  included  this  independence  does  not  occur,  as  will  be 
discussed  later;  however,  the  normalization  of  the  power 
density  is  used  in  all  regimes  for  continuity  and  for 
comparisons  with  Refs.  2-4. 

E.  EFFECTS  CP  HIND  VELOCITY 

As  wind  velocity  increases,  the  power  density 
distribution  spreads  out  and  the  maximum  power  density 
decreases  slightly.  This  effect  is  presented  in  Figs.  6-11 
for  various  angles  of  beam  incidence  with  respect  to  the 
vertical  (Psi  = 0,30,60  and  84  degrees).  Figures  10  and  11 
also  show  the  close  correlation  between  results  obtained 
from  the  Near  Zone  and  the  Far  Zone  solution  methods  for  Psi 
= 0 and  40  degrees,  respectively. 

The  spreading  effect  on  the  power  density  distribution 
is  due  to  larger  refraction  angles  (HR)  caused  by  higher 
wave  facet  slopes  that  occur  at  higher  wind  velocities.  The 
maximum  power  density  decrease  is  caused  by  the  spreading. 

Por  an  optical  communication  or  detection  system  an 
increase  in  wind  velocity  would  mean  that  the  maximum 
intensity  of  the  beam  that  could  be  focused  to  a point 
beneath  the  ocean  surface  would  be  reduced,  but  the  width  of 
receivable  signal  radiation  or  the  search  width  would  be 
increased. 

C.  EFFECTS  OP  INCIDENCE  ANGLE 


As  the  angle  of  incidence  with  respect  to  the  vertical 
(Psi)  is  increased  the  maximum  power  density  decreases  and 


the  distribution  shifts  away  from  the  vertical.  This  effect 
can  be  seen  in  Figs.  6-11  and  more  clearly  in  Figs.  12  and 
13.  Figure  12  is  a Far  Zone  solution  at  a depth  of  50 
meters  for  Psi  = 0,30,60  and  84  degrees  while  Fig.  13  is  a 
Near  Zone  solution  at  a depth  of  10  meters  for  the  same 
values  of  Psi. 

The  decrease  in  maximum  power  density  is  a result  of 
greater  reflectance  at  the  surface  facet  as  the  angles  of 
incidence  (HI)  increase  with  increases  in  Psi.  The  shift  in 
the  center  of  the  distributions  away  from  the  vertical  is 
caused  by  the  higher  angles  of  refraction  (HR)  that  occur  as 
HI  increases. 

For  an  optical  communication  or  detection  system  an 
increase  in  Psi  would  be  expected  to  lower  the  maximum  beam 
intensity  below  the  ocean  surface  and  to  shift  the  point  of 
maximum  intensity  away  from  the  point  of  water  entry. 


D.  EFFECTS  CF  HIND  DIRECTION 


The  effects  of  wind  direction  are  small  compared  to  the 
wind  velocity  effects.  Hhen  the  wind  direction  (Chi)  is 
perpendicular  to  the  beam  (Chi  = 90  degrees)  there  is 
slightly  less  spreading  of  the  power  density  distribution 
compared  to  Chi  = 0 or  180  degrees.  The  maximum  power 
density  is  essentially  unaffected  by  variations  in  Chi. 

Figure  14  shows  the  effect  of  variations  in  Chi  for 
normal  incidence  and  Fig.  15  for  oblique  incidence  (Psi  = 40 
degrees) . 

Hind  direction  appears  to  be  of  little  concern  to  an 
optical  communication  or  detection  system  as  related  to  the 


power  density  distribution  below  the  ocean  surface. 

E.  EFFECTS  CF  SPOT  SIZE 

Figures  16  and  17  show  the  effect  of  increasing  the 
incident  beam  spot  radius  from  0.1  to  0.5  meters  for  Psi  = 
0 and  40  degrees,  respectively.  This  is  done  with  the 
assumption  that  the  energy  per  unit  area  remains  constant; 
therefore,  the  total  energy  must  increase.  The  power 
density  spreads  out  and  increases  in  intensity  uniformly 
over  the  distribution  as  the  spot  size  is  increased. 

Since  the  equations  used  in  this  solution  method  deal 

with  the  ratio  of  P /P  , these  results  simply  mean  that 

d ds 

more  power  in  a larger  beam  will  produce  a more  intense  and 
more  widely  distributed  pattern  below  the  ocean  surface. 

F.  EFFECTS  OF  DEPTH 

The  decrease  of  power  density  with  increases  in  depth  is 
illustrated  in  Figs.  18  ^nd  19  for  Psi  equal  to  0 and  40 
degrees,  respectively.  The  decrease  is  linear  for  the 
normalized  power  density  as  can  be  seen  in  Fig.  20. 

Scattering  and  absorption  cause  this  decrease  in  the 
power  density.  As  the  depth  increases  the  diffuse 

transmittance  function  (DTF)  rapidly  becomes  an  inverse 
function  of  the  sum  of  the  hyperbolic  sine  of  the  depth  and 
the  hyperbolic  cosine  of  the  depth. 


▼ 1 

DTF  = f[  (sinhZ+coshZ)  ] 

which  means  the  DTF  decreases  rapidly  with  increases  in 
depth  and  approaches  zero  in  the  limit.  This  agrees  with 
the  theory  cf  Preisendorf er[ 8ef . 16]. 

Swennen's  prediction  of  a constant  normalized  power 
density  distribution  in  the  Far  Zone  regime  due  to  the 
diffusing  effects  of  the  ocean  in  Hefs.  2-4  holds  only  when 
the  effects  of  scattering  and  absorption  are  ignored.  One 
computer  run  was  made  with  the  DTF  set  egual  to  one  to 
illustrate  this  effect  in  Fig.  20. 

Figures  21  and  22  illustrate  the  decrease  in  power 
density  with  depth,  non-dimensionalized  with  respect  to  spot 
size,  Z/a,  again  for  Psi  egual  to  0 and  40  degrees, 
respectively . 

An  optical  communication  or  detection  system  will  be 
limited  greatly  by  the  power  density  decrease  associated 
with  depth. 


G.  MAXIMUM  POWER  DENSITY 


The  maximum  obtainable  power  density  available  over  the 
spectrum  of  incidence  angles  (Psi  =0-84  degrees)  for  a 
given  set  of  conditions  is  presented  in  Figs.  23  and  24. 
The  maximum  power  density  obtainable  does  not  occur  for 
normal  incidence,  but  for  a small  value  of  Psi,  because  of 
the  fact  that  as  the  wind  increases  the  most  probable  slope 
is  not  zero  but  a small  angle.  The  maximum  obtainable  power 
density  then  drops  continuously  with  increases  in  Psi  after 
the  peak  that  occurs  around  Psi  =10-15  degrees. 
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Figure  23  illustrates  the  decrease  in  the  naximum 
obtainable  power  density  for  all  values  of  Psi  as  the  wind 
velocity  increases.  This  is  due  to  the  spreading  effect  on 
the  power  density  distribution  caused  by  the  wind  that  was 
described  previously. 

Figure  24  illustrates  the  independence  of  the  naximua 
power  density  on  the  wind  direction  (Chi)  for  all  values  of 
Psi.  Chi  egual  to  270  degrees  is  not  plotted  but  is 
equivalent  to  the  condition  at  90  degrees  because  of  the 
symmetry  of  the  slope  probability  function  about  the  wind 
direction. 

The  angle  with  respect  to  the  vertical  from  the  point  of 

observation  below  the  ocean  surface  (Hu  ) at  which  the 

0 

maximum  power  density  occurs  for  each  value  of  Psi  is 

plotted  in  Fig.  25.  The  plot  is  independent  of  wind 
velocity  and  direction.  Examining  Figs.  16  - 19  it  can  be 
seen  that  Fig.  25  is  also  independent  of  depth  and  spot 
size.  Mu^  for  maximum  power  density  increases  nearly 

linearly  from  approximately  1 to  56  degrees  as  Psi  goes  from 

0 to  84  degrees.  The  increase  in  Hu^  is  due  to  the  power 

density  distribution  shifting  away  from  the  vertical  for 

increasing  valules  of  Psi  as  discussed  earlier  and 
illustrated  in  Figs.  12  and  13. 

The  obvious  consideration  with  an  optical  communication 

or  detection  system  is  the  decrease  in  maximum  power  density 

that  accompanies  higher  incidence  angles  and  wind 

velocities.  The  strong  independence  of  the  angle  Mu  at 

0 

which  the  maximum  power  density  occurs  with  respect  to  the 
wind  velocxty,  wind  direction,  depth  and  spot  size  indicates 
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VI . CONCLUSIONS 


The  time-averaged  power-density  distribution  below  the 
ocean  surface  due  to  incident  laser  radiation  has  been 
examined  through  computer  simulations  in  both  the  Near  Zone 
and  the  Far  Zone  regimes.  The  effects  of  altering  incidence 
angles,  wind  velocities,  wind  directions,  beam  spot  sizes 
and  depths  of  penetration  have  been  presented. 

The  power-density  distribution  was  found  to  be  highly 
dependent  on  the  angle  of  beam  incidence,  the  wind  velocity 
and  the  attenuation  associated  with  depth  of  penetration  of 
the  laser  beam.  The  fact  that  the  location  of  the  maximum 
power  density  was  found  to  be -fundamentally  dependent  on 
only  the  angle  of  beam  incidence  indicates  a possibility  of 
direction  finding  capabilities  by  a submerged  receiver. 

Improvements  in  computer  systems  and  analysis  techniques 
have  allowed  an  examination  of  a much  broader  range  of 
variables  than  were  possible  heretofore.  Angles  of  beam 
incidence  up  to  84  degrees  were  considered  within  the  range 
of  validity  fcr  restrictions  placed  on  the  simulation. 
Above  84  degrees  incidence  angle  the  possible  interactions 
at  the  air/sea  interface  make  an  analysis  of  the 
power-density  distribution  extremely  complex,  and 
experimentation  will  be  necessary  to  give  indications  of  the 
feasibility  of  laser  operations  in  this  area. 

This  simulation  could  be  extended  to  an  analysis  of  much 
larger  spot  sizes  that  would  be  expected  from  a satellite 
based  transmitter.  Prettyman  and  Cermak[Ref.  17]  suggested 
that  the  effects  of  larger  spots  could  be  approximated  by 


assuming  each  wave  facet  independent  of  its  neighbors  and 
then  integrating  over  the  larger  spot  area,  since  the  small 
spot  represents  the  outer  limit  of  the  effects  of  the  ocean 
surface  on  a laser  beam  of  greater  spot  size.  There  would 
have  to  be  a large  number  of  surface  facets  contained  in  the 
spot  to  insure  facet  independence  and  accurate 
approximations,  hence  a strong  dependence  on  the  wavelength 
of  the  wave  associated  with  wind  velocity. 

An  extension  of  the  simulation  to  pulsed  laser  radiation 
is  also  warrented,  since  this  mode  of  operation  may  be 
required  for  air-to-subsurface  laser  systems.  The  effects 
of  pulsed  radiation  occur  nearly  instantaneously  as  compared 
to  the  time-averaged  analysis  performed  in  this  simulation. 

A further  refinement  of  the  simulation  would  involve  a 
more  complex  modeling  of  the  beam  attenuation  due  to 
scattering  and  absorption  in  the  ocean  through  the  use  of  a 
diffusion  or  a multiple  forward  scattering  method,  both 
described  by  BucherfHef.  18]  and  Preisendorf  er[  Ref . 16]. 
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APPENDIX  B 


COMPOTES  PSOGBAM  OSAGE 


To  use  the  computer  programs  listed  following  this 
Appendix,  first  compute  the  ratio  of  the  depth  of  the  point 
of  observation  below  the  ocean  surface  to  the  beam  radius  at 
the  surface  facet,  Z/a.  Ose  the  Near  Zone  program  for 
Z/a<100  and  the  Far  Zone  program  for  Z/a>100.  The  Near  Zone 
program  may  be  used  for  any  value  of  Z/a,  but  a slight 
computation  time  and  computer  storage  increase  will  result. 

The  only  other  basic  reguirement  is  the  selection  of 
eight  input  variables.  The  variables  are  listed  below  along 
with  the  line  from  the  respective  program  in  which  they 
appear. 


NEAR 

FAR 

VARIABLE 

ZONE 

ZONE 

Depth,  Z (m) 

470 

44  0 

Beam  radius,  RAD  (m) 

480 

450 

Angle  Nuq,  NUORON  (deg) 

490 

480 

Backward  scattering  coefficient,  BS 

560 

530 

-1 

Absorption  coefficient,  AC  (m  ) 

570 

540 

Beam  incidence  angle,  PSIR0N  (deg) 

630 

590 

Hind  direction,  CHI  (rad) 

680 

640 

Hind  velocity,  W (m/sec) 

710 

670 

Z,  RAD,  NOORUN,  BS  and  AC  are  direct 

program 

inputs 
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aade  on  the  indicated  lines.  PSIr  CHI  and  H are  selected 
from  DATA  statements  430  and  440  in  the  Near  Zone  program 
and  420  and  430  in  the  Far  Zone  program.  The  variable 
locations  in  DATA  statements  mere  used  for  multiple 
variations  of  the  parameters  in  a slightly  different  form  of 
the  program  than  presented  which  is  easily  set  up  with  DO 
LOOPS.  The  entry  variables  may  also  be  varied  in  this  same 
manner. 

Both  programs  must  be  combined  with  the  common  package 
of  subroutines  that  are  listed  following  the  basic  program 
listings. 

The  output  from  both  listed  programs  appear  on  the  next 
two  pages. 
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